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ABSTRACT 

Aims. We present a detailed X-ray and radio wavelength study of G308.4-1.4, a candidate supernova remnant (SNR) in the ROSAT 
All Sky Survey and the MOST supernova remnant catalogue. 

Methods. The SNR candidate and its central point sources were studied using observations from the Chandra X-ray Observatory, 
SWIFT, the Australian Telescope Compact Array (ATCA) at 1.4 and 2.5 GHz and WISE infrared observation at 24 /jm. 
Results. We conclude that G308.4-1.4 is indeed a supernova remnant by means of its morphology matching at X-ray, radio and 
infrared wavelength, its spectral energy distribution in the X-ray band and its emission characteristics in the radio band. G308.4-1.4 is 
a shell-type SNR. X-ray, radio and infrared emission is seen only in the eastern part of the remnant due to a strong spatial density 
variation of the interstellar medium around the remnant. The X-ray emission can best be described by an absorbed non-equilibrium 
collisional plasma with a hydrogen density of = (1.02 ± 0.04) cm"^, a plasma temperature of 6.3*^ ^ million Kelvin and an under- 
abundance of Iron, Neon and Magnesium, as well as an overabundance in Sulfur with respect to the solar values. A Sedov analysis 
revealed that the remnant is at a distance of « 10 kpc and the progenitor star exploded ~ 5000 to 7500 years ago. Two faint X-ray 
point sources located near to the remnant's geometrical center are detected. Both sources have no counterpart at other wavelengths, 
leaving them as candidates for the compact remnant of the supernova explosion. 

Key words. ISM: supernova remnants - ISM: individual objects: G308.4-1.4- Stars: neutron 



■ 1. Introduction 

' Supernovae (SNe), the core-collapse of massive stars or the ther- 
. monuclear disruption of white dwarfs, are r are events, occur- 

■ ring o n average every 50 years in the Galaxy (iKeane & Kramer! 

■ 12008b . However, in the past 2000 years, only 7 Galactic SNe 
I have been observed - SN 185 (RCW 86), SN 386 (Gl 1.2-0.3), 

■ SN 1006, SN 1181 (3C58), Crab SN, Tycho SN and the Kepler 
[ SN. Most Galactic SNs are unobserved owing to visible-band 

extinction by interstellar dust. Mor eover, only ^ 274 SNR have 
been identified in the radio band (i GreenI 120091 and references 
, therein). Although the Green catalog represents the results of 
more than 40 years of intensive search with the largest radio 
i telescopes, it is incomplete and strongly biased by selection ef- 
. fects. For a source to be included in a supernova remnant cata- 
' log, it must be (1) extended, with a diameter > 3'; (2) located 
I at low galactic latitude; and (3) have a non-thermal spectrum 
. with a spectral inde x in the centimeter-band steeper than 0.3 
; (lHelfand etal.1119891). 

Most SNRs emit thermal radiation when the SN blast wave 

■ expands in the surrounding interstellar medium (ISM), form a 
shock wave at the shock front which then ionizes the atoms 
and increases the temperature to 10^ - 10^ K. Therefore, the 
successful completion of the ROSAT All-Sky Survey (RASS) 
in the X-ray regime provided a new and exciting window 
for both finding SNRs and the compact objects that may re- 
side within them. The seven X-ray detected and radio-quiet 
Central Compact Objects which were associated with SNRs (see 
iGotthelf & Halpernll2008l for a review) are a good example for 
the impact on pulsar and SN science of X-ray observatories. 
Additionally, X-ray observations of SNRs at moderate spectral 
resolution can determine remnants properties like shock veloci- 
ties, post-shock gas densities and temperature as well as swept- 
up shell mass and overall morphology. 



iBussed (Il998h and ISchaudell (l2003h have demonstrated 
the potential power of X-ray surveys for finding new SNR 
candidates. They analyzed ROSAT all-sky survey data in 
a systematic search for extended (> 5') X-ray sources at 
low Galactic latitude (b < +15°) by correlating them with 
databases like SIMBAD, NED, NVSS, SUMSS, NRAO, ATNF, 
Parkes and Effelsberg radio survey data. Finally, they cor- 
related the digital sky survey (DSS II) images for bright 
optical sources close to the candidate's positions. Of their 
candidate list 14 sources have been identified meanwhile 
as SNR (Pfeffermann& Aschenbach 1996; 'Reich et al." 19921 
Busser et al. 1996; Asaoka & Aschenbach 1994; AsaokaetaL 
1996; Eggeretal. 1996; Folaheraiter et al. 1996; Jackson et al. 
2008: .Robbins etal.. ,2011; Stupar et al..,2007i: iTian et al. 201oT 
These identifications suggest that many of their SNR candidates 
are indeed radio-under-luminous, explaining why past radio sur- 
veys could not detect or identify them before. 

In a further down-selection of promising candidates our list 
contains now 9 sources which appear to be very promising SNR 
ca ndidates. One of them is G308.4-1 . 4, list ed as G308.3-1.4 
in lSchaudeleTaP (l2002h and lSchaudei (l2003h . Since the RASS 
had only an average exposure time of ~ 400s for sources in the 
Galactic plane the available ROSAT data do not support a de- 
tailed spectral and spatial analysis of these sources. However, 
G308.4-1.4 appears in the RASS data as a spherical center-filled 
X-ray source whereas in the Sydney University Molonglo Sky 
Survey (SUMSS) two radio arcs have been detected which partly 
overlap with the RASS data (see figure [B. Moreover, G308.4- 
1 .4 is listed in the MO ST supernovae remnant catalogue as pos- 
sible SNR candidate (.Whiteoak & Greenll 19961) . 

In this paper we report on the identification of the SNR 
candidate G308.4-1.4 as SNR using archival Chandra, SWIFT, 
ATCA and WISE observations. The results of the spatial and 
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Fig.l. ROSAT RASS 30' x 30' image of G308.4-1.4 (0.1-2.4 
keV) overlaid with the MOST survey image at 843 MHz (con- 
tour lines at 4, 8, 12, 16 and 25 mJy). The half-power beam- 
width of the telescope is 43". The RASS image is binned with 
0.75' per pixel and smoothed by a Gaussian kernel of cr a; 2'. 
The color bar at the right specifies the color assignment in units 
of counts/image pixel. 



spectral analysis of the X-ray and radio data are presented in ^ 
and |3j respectively. In ^we report on data of G308.4-1.4 in 
other wavelength regimes. A summary and a concluding discus- 
sion about the identification of G308.4-1.4 and its central point 
sources is given in ^ Furthermore, in the same section we use 
the inferred spectral parameters of the SNR to derive fundamen- 
tal parameters of the remnant. 



2. X-ray observations and data analysis 

The SNR-candidate G308.4-1.4 was observed with the 
Advanced CCD Imaging Spectrometer (ACIS-I) on board the 
Chandra X-ray observatory on 26th and 27th of June 2010 
(obs.ID. 11249) for x 15.1 ks. All data were uncontaminated 
by soft proton flares. After dead-time correction we were left 
with 14.9 ks of good data. Standard processed level-2 data were 
used with the newest de-gap and tap-ringing corrections applied. 
The reduction of the data was done with the Chandra Interactive 
Analysis software CIAO, version 4.3. For the spectral fitting we 
used the XSPEC software package, version 12.7.0. 

In addition to Chandra we made use of the SWIFT observa- 
tory to analyze X-ray point sources detected in the central re- 
gion of the remnant. G308.4-1.4 was observed three times with 
the Swift X-ray Telescope (XRT) on the 17th of June (obs.ID 
00032030001), the 7th (obs.ID 00032030002) and the 11th of 
August (obs.ID 00032030003) for 3.5 ks, 4.2 ks and 1.3 ks, re- 
spectively. These data were reduced with the HEASOFT soft- 
ware package, version 6.1 1. For the spatial and spectral analysis 
we use the high resolution Chandra data as it contains five times 
more counts of G308.4-1.4. 

To calculate the error in the X-ray flux we used the XSPEC 
tool fluxerror.tcQ. For the analysis of the extended source emis- 
sion we restricted the energy range to 0.5 - 4.0 keV. All errors 
in this publication represent the Icr confidence range for one pa- 
rameter of interest, unless stated otherwise. 
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Fig.2. Chandra ACIS-I color image of G308.4-1.4 (red 0.1- 
1.0 keV, green 1-1.5 keV and blue 1.5-5 keV). The images are 
binned with 1" per pixel and smoothed by a Gaussian kernel of 
cr = 3". 20 point sources are detected in the 12' x 12' area around 
the center of G308.4-1.4 (c.f. table©. 



2.1. Spatial analysis 

Inspecting Fig.l, the X-ray image of G308.4-1.4 shows strong 
emission in the eastern part. This region has a circular shape with 
a diameter of ~ 8.4 arcmin (see fig.|2]i. The radiation is fading 
towards the geometrical center of the source. The coordinate of 
the circle center is RA,. 13Hl'"24.9^ DEC, ~ -63°43'07.8" 
which was obtained by fitting an annulus to the outer skirts of 
the extended emission. 

Searching for point sources in a 12' x 12' box around 
G308.4-1.4 using the wavelet source detection algorithm 
wavdetect we found 21 sources. Their position, positional er- 
ror and signal-to-noise ratio is listed in table [T| and marked by a 
circle in figure |2] Cross-correlating the sources with the help of 
the VizieR online tool we found seven matchesQ These are listed 
with the positional discrepancy in table [T] 

Eight of the detected sources are located within the ex- 
tended emission of G308.4-1.4 and two of them (# 1 and # 10) 
are within ^ 1 arcmin of the geometrical center of the SNR- 
candidate. Source 1 is the brightest among all detected point 
sources and its designation is CXOU J134124.22-634352.0. 
However, its position is coincident with 2MASS J13412422- 
6343520. In the optical source catalogs this source has no in- 
formation about its proper motion. Hence we cannot put better 
constraints on the positional coincidence of source # 1 and the 
2MASS object. Source 10 has no optical counterpart and is des- 
ignated with CXOU J 1 34 1 27 . 1 2-634327.7. 

2.2. Spectral analysis of G308.4-1.4 

To extract the spectra of the SNR-candidate we first removed 
all point like sources and then selected all photons in a circular 
sector with a radius of 4.2 arcmin and central angle of ^ 136° 



' |http : / /heasarc . nasa . gov/xanadu/ xspec/f luxerror . html | ^ |http : //vizier .u-strasbg. fr/viz-bin/VizieR| 
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Table 1. Detected sources in a 12' x 12' box around G308.4-1.4. The sources are denoted as in figure|2l 



Source 


RA (J2000) 


DEC (J2000) 


SRA 


<5DEC 


S/N" 


optical counterpart (positional discrepancy) 




h:m:s 


d:m:s 


arc sec 


arc sec 






1 


13:41:24.217 


-63:43:52.04 


0.41 


0.40 


63.0 


2MASS J13412422-6343520 (0.07') 


2 


13:40:50.534 


-63:38:03.12 


0.45 


0.42 


2.7 




3 


13:41:24.203 


-63:47:24.47 


0.70 


0.42 


3.1 




4 


13:41:21.634 


-63:47:05.63 


0.44 


0.43 


25.5 


2MASS J 13412162-634705 3 (0.27') 


5 


13:40:38.921 


-63:46:38.83 


0.54 


0.45 


3.3 


2MASS 13403887-6346384 (0.48') 


6 


13:41:05.585 


-63:46:25.89 


0.47 


0.42 


4.9 




7 


13:41:10.742 


-63:45:36.11 


0.52 


0.42 


4.9 




8 


13:41:04.580 


-63:44:56.66 


0.52 


0.42 


6.7 




9 


13:40:44.230 


-63:44:29.50 


0.65 


0.41 


2.4 




10 


13:41:27.121 


-63:43:27.72 


0.56 


0.44 


3.8 




11 


13:41:14.267 


-63:39:45.91 


0.57 


0.43 


9.7 


3UC 053-222814 (0.02') 


12 


13:40:38.666 


-63:39:04.77 


0.61 


0.43 


9.0 


2MASS J13403862-6339044 (0.41') 


13 


13:41:02.437 


-63:38:00.00 


0.56 


0.44 


0.8 




14 


13:40:36.217 


-63:45:57.41 


0.47 


0.44 


2.0 




15 


13:41:01.308 


-63:43:23.50 


0.40 


0.40 


2.9 




16 


13:41:11.249 


-63:39:44.07 


0.76 


0.44 


4.7 




17 


13:41:54.767 


-63:39:20.68 


0.91 


0.67 


6.7 


2MASS J13415473-6339207 (0.22') 


18 


13:41:18.147 


-63:37:01.87 


1.31 


0.53 


3.8 




19 


13:41:14.934 


-63:48:38.00 


0.67 


0.54 


3.2 




20 


13:41:49.497 


-63:46:23.59 


1.06 


0.65 


4.1 




21 


13:41:43.058 


-63:38:38.10 


1.24 


0.50 


5.4 


2MASS J13414304-6338371 (0.93') 



" ctg = 1 + -^c'bg -I- 0.75, where Cbg are the background counts. 



(see figure |2]). The background spectrum was derived from all 
photons in an annular sector with the same central angle as the 
circular one but with the radii 5.5' and 7.2'. After subtraction 
of the background contribution (0.054 cts/pix) we were left with 
1 1 752 counts for the spectral analysis of the whole remnant- 
candidate. Additionally, the emission region of G308.4-1.4 was 
divided into three separate regions which are illustrated in figure 
|2l The net counts are 4095 for the inner sector and 7658 for the 
outer shell, as well as 4499 for the northern part of the outer shell 
and 3158 for the southern part. The spectrum of every single re- 
gion was binned with at least 50 counts per bin. Only the spectra 
of the southern and northern parts of the outer rim were binned 
with at least 75 counts per bin. 

To fit the X-ray spectrum of G308.4-1.4 we tried var- 
ious one- and two-component spectral models. However, 
only a model of an absorbed collisional plasma which 
is in non-equilib rium and allows a temperature evolution 
dBorkowski et alJ2 001) reproduces the energy distribution of the 
extended source. In XSPEC this model is implemented as vgnei. 
To improve the goodness of the fit and to obtain the abundance 
of metals in the plasma we thawed the parameter of every single 
metal and monitored the improvement of the fit statistic. 

The analysis results in a best-fit model with a reduced 
of 1.03 for 95 degrees of freedom (d.o.f.) and a hydrogen col- 
umn density hh - 1-02 + 0.04 x 10^^ cm"^. The value for hh 
is smaller than the integrated hydroge n column densi ty towards 
G308.4-1.4whichisL31><10f_cm-2 (Kalberla et al. 2005) and 
1.45 X 10^^ cm"^ (iDickev & Lockman 1990) . The temperature 
of the plasma is found to be 6.3;*^g y million Kevin. The best-fit 
was obtained when thawing the abundance in Neon, Magnesium, 
Sulfur and Iron which differs significantly from the solar values: 
Ne = (0.72 + 0.14) Neo, Mg = 0.70!^;^^ Mg^, S = 2.3+°;^ So and 
Fe - 0.70;^q [2 Fe©. The ionization timescale tq = foA^e ik) is the 
remnant's age and A^e the post-shock electron number density) is 
2.7 X 10'" s/cm-' and the ionization timescale-averaged tempera- 
ture {kT) = f'° T(t)N,(t)dt/T = 1.2+2-2 j^gY jjjg jjjj^g jjjg 



shock). Thus, the un-absorbed flux in the energy band 0.5 - 4.0 
keV of G308.4-l.4is fx = 1.61+°-^^ x lO"'" ergs cm^^s-i. The 
model fit and the raw spectrum is shown in figure[3]and the fitted 
parameters are listed in table|2] 
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Fig. 3. Spectrum and fitted model of the X-ray emission of the 
whole remnant. 



Adding a second gnei component results \n x^^^ - 106 for 
91 d.o.f. and hence is no improvement over a single model. The 
other models we have tried to fit to the spectral distribution are 
a VMEKAL-model (xl^^ = 2.92 for 96 d.o.f.), a VAPEC-model 
{X^.^ii - 2.34 for 96 d.o.f.) and a thermal bremsstrahlung model 
^rerf ~ ^-^ ^'^^ d.o.f.). Non of these models can reproduce 
the spectrum of G308.4-1.4. However, all three model fits re- 
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suits in a hh ~ 1 X 10^^ cm ^ and a plasma temperature of ~ 0.5 
keV. 



2.2.1 . The inner sector and the outer rim 

To investigate the spatial variation of the hot plasma we have 
divided the extended source into two regions, the inner sector 
and the outer rim. Additionally, the outer rim was divided into 
two parts, north and south (see figure |2] for a illustration). For 
the background spectra of the four regions we used the same 
area on the detector as for the whole remnant, only for the outer 
rim north and south we used a center angle that matched the one 
from the source regions. 

First we fitted the plasma model with free ria. In every of the 
four spectral fits we found a hydrogen column density which is 
in good agreement with the rtu of the whole remnant. Second 
we fitted the spectra with a ne fixed to the value of the whole 
remnant. The fitted value for tin and the results of the fits with 
a frozen riu are summarized in table |2] and the spectra with the 
respective fits are shown in figureH] 

Within the statistical uncertainties, all examined regions turn 
out to have a similar temperature. Only the ionization timescale- 
averaged temperature is varying. The large statistical uncertainty 
also prevents us from concluding on a temperature evolution 
with the distance to the expected center of the remnant. No vari- 
ation of the Neon and Magnesium abundance over the surface 
of the remnant is detectable too. The only noticeable change 
in abundance is in Sulfur and Iron in the outer rim. The spec- 
trum in the northern rim shows an abundance in Iron compa- 
rable with the solar value. In the southern rim no significant 
deviation from the solar value was found in the abundance of 
Sulfur The ionization timescale varies with distance from the 
remnant's center and the direction of the expanding shock wave 
into the ISM. The highest values for r were found for the in- 
ner sector (8.1 x 10'° s/cm^) and the lower one in the southern 
part of the rim (2.0 x 10'" s/cnr') which implies a different post- 
shock electron number density in the inner and outer regions of 
the SNR-candidate. 



2.3. The central point sources 

As mentioned above, two point sources were detected 
close to the presumed center of the remnant, source # 1 
(CXOU J134124.22-634352.0) and # 10 (CXOU J134127.12- 
634327.7). For the spectral analysis of source # 1 we used all 
photons within a circular region of radius 10 arcsec centered on 
the source position. The background was extracted from an an- 
nulus with radii 15" and 20". After background correction we 
were left with 195 source counts in the energy band between 
0.5 and 10.0 keV. We used the tool psextract to compute the re- 
sponse matrices and binned the spectrum with at least 20 counts 
per bin. 

A single absorbed blackbody or a single absorbed power law 
model for fitting the spectral distribution gives poor results. In 
the case of a blackbody spectrum the reduced - 2.67 for 6 
degrees of freedom is unacceptably large. The power law model 
fit results in a reduced of 0.94 (6 d.o.f.). However, the post- 
fit residuals show systematic derivations between the data points 
and the fitted model. Therefore, we used a two-component ab- 
sorbed blackbody model which results in - 1-62 (4 d.o.f.) 
and a hh = 1.34;^[ 3° x 10^^ cm"^. In order to better constrain the 
spectral parameters we fixed ne to the hydrogen column density 
obtained from fits to G308.4-1.4. This gave a reasonable fit with 




Distance [kpc] 

Fig. 5. Dependence of the blackbody emitting radius Remii on the 
distance to source # 1 in kpc. The error range is gray shaded and 
the green fine indicates the upper limit on the radius of a neutron 
star. 



- I -96 for 5 d.o.f. (see figure|6]l. The temperatures of the two 
blackbodies are Ti = 1.3!°-^ x lO*" K and T2 = 6.4+2-^ x 10^ 
K. Using the normalization factors Ki = 7;^^' x 10"^ and K2 = 
(1.3 + H-0.4) X 10"^ of the two blackbodies and the distance to 
the source Dio in units of 10 kpc we estimate the emitting source 
radius R as follows: 



/Txio^y ^ ^ 
^ - 1 1 — 7 M, 



(1) 



where ctb is the Stefan-Boltzmann constant. As we do not know 
the distance to G308.4-1.4 we consider the emitting radius as 
function of the distance to the source. This is illustrated in figure 
|5] The result for the second blackbody is R2 = 340;^jg^^ii^J m. 
The un-absorbed flux in the energy band 0.5 to 10 keV is fx = 
2.2^j g X 10"'^ ergs cm^^s"' and in the energy band 0.3 to 3.5 
/x = (4 + 2) X 10"'^ ergs cm^^^-i 

Additionally, we tried a model composed of a blackbody 
and a power law with the same fixed «h as before. This gives 
a reasonable fit with;i'^ = 2.17 (5 d.o.f.). The photon index is 
2.4+[ '2' and the blackbody temperature is (1.3 + 0.3) x 10^ K. 
With the normalization of Ki = 9^^^ x 10"^ the emitting radius 
is 75^4g''t/g g km. The un-absorbed flux in the energy band 0.5 to 
10 keV is fx = 2.5^J^l x 10"'^ ergs cm^^g-i. 

As mentioned in section 12.11 the position of source # 1 is 
consistent with 2MASS J 134 12422-6343520, a star with a mag- 
nitude of B = 19.7, R = 19.6, J = 14.0, H = 13.3 and K = 
13.1. This implies an X-ray-t o- visual flux ratio of lo g( fx I fy) = 
log(/x) + V/2.5-H5.37 = 1.55 (iMaccacaro et al.ll988h . Here fx is 
the X-ray flux in the 0.3-3.5 keV band in ergs cm"'^ s ', V is the 
visual magnitude and we assumed that V ^ R. Fo r stars the ratio 
shoul d be in the range \og{ fxlfv) - -2.46 + 1.27 dKrautter et al.l 
Il999h . 

In addition we tried to fit various plasma models to the 
source spectrum. The only reasonable fit was found when using 
a APEC model. The goodness of fit which we found for 6 d.o.f. 
was x]-ed ~ '^■^^ ^'^^ '-^^ following spectral parameters: hh < 
0.47 x\(fi^ cm-2, T = 3.1!|-^ x 10^ K. For this model the flux 
in the 0.3 to 3.5 keV range is fx - l-l^oi x 10'^ ergs cm^^s"'. 
Thus, the X-ray-to-visual flux ratio is \og{fxlfv) - 0.25, still 
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Table 2. Detected sources in a 12' x 12' box aroutid G308.4-1.4. The sources are denoted as in figure|2l 





Whole remnant 


Inner rim 


Outer rim 


Outer rim north 


Outer rim south 


nu riO-^cm"-l 


1.02 ± 0.04 


1 05+"" 


1 00 + 05 


1 02 + 09 


0.99 + 0.9 


k^T [keV] 


0.54+«'« 

-0.06 


0.55+"-?^ 

-0.10 


0.54+""? 

-0.05 


0.50+"'2 

-0.06 


0.57 ± 0.08 




0.72 ± 0.14 


0.7 + 0.3 


0.73+"!^ 

" -^-0.14 


0.5 ± 0.2 


0.8 ±0.2 




0.70+""? 

-0.08 


0.80+"?? 

" " -0.17 


0.67+"'" 

" " -0.09 


0.83+"ll 


0.62+"!i 

-0.12 




9 Q+0.7 
^■-'-0.8 


3 0+' ' 

-'•"-1.0 


1 7+" 


2 4+1.3 




Fe [Feo] 


70+° '" 

' -0.12 


0.90 ± 0.3 


64+" '" 

"■"^-0.12 




45+"" 


T [10'" s/cm^] 


2.7 


8.1 


2.1 


5.2 


2.0 


{kT) [keV] 


1 2+2-2 

^■^-0.4 


72+" " 

"■'^-0.12 


1 4+1.6 




1 1+'" 


Norm [10"^ cm"']" 


3.3 ± 0.9 


9+"-^ 

"■^-0.2 


2.4 ± 0.5 


1.4 ±0.4 


1.0 ±0.3 


^* 


98.19 


57.54 


76.81 


38.74 


18.55 


d.o.f. 


95 


57 


77 


38 


27 


/x (0.5-4.0 keV) [lO"'" ergs cm-^s-'] 


1 61+" '^ 


30+" "' 

"•■^ -0.02 


1 9^+0.08 
-^-0.09 


55+" "2 

"•-^-^-0.04 


56+" "^ 

"••^"-0.06 



" Norm= ,n"n"' >i2 f NcNndY, where Da is the angular diameter distance to the source in cm, A'e and Nn are the post-shock electron and hydrogen 
densities in cm"^, respectively. 




Energy (keV) Energy (keV) 



Fig. 4. Spectrum and fitted model of the X-ray emission of the various regions. Upper left: inner part of the remnant, upper right: 
outer rim, lower left: outer rim north and lower right: outer rim south. 



higher than the normal value for stars. Again no variation in the are detectable. All fitted models and their inferred spectral pa- 
spectral parameters between the flaring and quiescence epoch rameters are listed in table [3] 
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Table 3. Spectral fit results of the central point source # 1 . 



Model 


X^I±o.i. 


"H 

10^2 cm-' 


TIkT, 
-/keV 


Ki 


kT2 

keV 




f 0.5-10/ f 0.3-3.5 

Jx iJx 
10"'2ergs cm-2 s"' 


PL"- 
BB+BB^ 
PL+BBA 
APEC 
BREMS'"' 
RAYMOND" 


16.00/6 
5.62/6 
1.96/5 
2.74/5 
4.66/6 
6.98/6 
5.70/6 


1.2 X 10-" 
<0.33 

1 ^4+1-80 

< 3.8 

< 0.47 
1.8 X 10-' 

< 0.21 


0.44 

9 9+0.9 
^•^-0.4 

9 4+1.0 
2 7+'-5 

^•'-1.2 

2 1+2-1 
2-8!u 


9.9 X 10-^ 
3.9+_ll X 10-5 

7+21 X 10-5 

7+^1 X 10-5 
9.0![§ X 10-5 
5. 1^2 3 X 10-5 
8.9!i;^ X 10-5 


55+''-24 
0.11 ±0.03 


(1.3 + -0.4) X 10-^ 
9+f X 10-5 


-/- 

16+'"'V0 16+"'"' 

"•^"-0.02'"-^"-0.03 

2.2+^-5/4 ± 2 
2.5!«;5/5!3 

^ -'-0.02' "-^ ^-0.03 
-/- 

0.13±-0.02/0.11!°;|!2 



" fit shows systematic derivations between the data points and the fitted model. 

the spectral parameters are calculated with nn fixed to the value obtained in fits to G308.4-1.4. 

Power law norm in photons keV-lcm-2s-' at 1 keV. 
" some errors could not be calculated. 




1 - 







-1 



1 2 
Energy (keV) 

Fig. 6. Spectrum of source # 1 fitted with a double blackbody 
spectrum. 



2.3.1 . Temporal analysis of source # 1 

In order to search for any temporal variability we binned the 
photon time of arrivals with 400 s per bin (see figure |7]i. The 
data show clearly a flare in the first ks of the observation. The 
spectral parameter of the double blackbody in the flaring and 
quiescence epoch match the overall spectrum within the error 
bars. Varying only the normalizations of the two blackbodies the 
flux in the quiescence epoch is with /_^*'^-'" = (1.9 + 0.3) x 10-'^ 
ergs cm~^s~' slightly lower than the flux in the total exposure. 
However, the flux in the time interval when the source was flar- 
ing is /;^'''-'° - (5.2 + 1.5) X 10-'^ ergs cm-^s-', three times 
the quiescence flux. However, the sparse photon statistics lead 
to large uncertainties in the derived errors of the quiescence and 
flaring fluxes and thus they are only rough estimates. Merging 
the SWIFT observations we cannot detect a source at the posi- 
tion of source # 1 . In the merged image a total of 4 counts were 
recorded within a circle with radius 15 arcsec (equals an encir- 
cled energy of 70 %) centered on the position of source # 1 . The 
total exposure time of the source was 8344 s. Hence, the 3cr up- 
per limit for the counting rate is approximately 8.6 x 10-"^ cts/s. 
Using WebPIMMS with the SWIFT response and the spectral 




, , , , \ , , , , [1 , , , , \ ^ i- , , 

55374.00 55374.05 55374.10 55374.1^ 

Time [MJD] 

Fig. 7. Flux variation of source # 1 in the Chandra observation 
binned with 400 s per bin and scaled with the mean flux in the 
whole exposure. The red arrow indicates the 3cr upper limit de- 
rived with the SWIFT observations. 



parameters derived for a double blackbody spectrum the upper 
limit on the flux is ^ 5.3 x 10-'^ ergs cm-^ s-'. For the APEC 
model the upper limit is ^ 2.3 x 10-'^ ergs cm-^ s-' . These limits 
are less than 25 % of the average flux in the Chandra observa- 
tion, clearly indicating that the source flux is variable on time 
scales of at least hours to weeks. 

2.3.2. Extent of source # 1 

With the help of a one-dimensional distribution of the 
counts/0.5" along a vertical line in the Chandra image around 
source # 1 we checked if the source is extended. The resulting 
diagram was obtained by integrating all counts within 5" x 0.5" 
rectangular apertures moving along the western direction with 
0.5" steps. We compared this with the count distribution of a 
point source simulated with MARX and the spectral informa- 
tion above0 However, no excess emission of point source # 1 
could be detected. 



^ See |http : // space .mit . edu/ cxc/marxT] 
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For source #10, the other source located close to the center 
of the remnant, we detected only 8 net counts. Hence, the count 
rate is (5.1 + 1.9) x 10""* cts/s, but no spectral analysis is possible 
by the Umited statistics. 

3. Radio observation and data analysis 

The source G308.4-1.4 was observed for 11.25 hours with the 
Australia Telescope Compact Array, a synthesis telescope near 
Nan-abri, New South Wales on 11th January 2002 at 1.432 and 
2.448 GHz with 13 channels per band. The single pointing con- 
tinuum observations of the source G308.4-1.4 was carried out 
with the array configuration 0.75C (maximal baseline length is 
750 m). All the four Stokes parameters have been recorded. 

The flux density calibration was performed through observa- 
tions of PKS B1934-638, which is the standard primary calibra- 
tor for ATCA observations. Phases were calibrated using obser- 
vations of secondary calibrating sources PKS 1329-665. Since 
the primary beam response is frequency dependent we did not 
merge the data from two observing bands before imaging and 
cleaning. Every single observing run and each of the two ob- 
serving bands were calibrated separately following the standard 
procedures for ATCA observations. 

The reduction was carried out using the ATNF release of 
the Multichannel I mage Reconstruc tion, Image Analysis and 
Display (MIRIAD. ISault et al.lll995h . A number of steps in the 
reduction process (i.e. "flagging") were done interactively. 

Because of the bad phase stability during the observation, 
the scalar averaging of gains over the interval of 5 min (the 
length of secondary calibrator's scan) was performed. Then, a 
total-intensity continuum image square region was formed using 
multi-frequency synthesis, uniform weighting and a cell size of 
10 arcsec. At the stage of the dirty map production all the corre- 
lations with 6th baseline were excluded to obtain an image with 
better signal-to-noise ratio. The ima ge was then de-convolved 
using the standard CLEAN algorithm (Clark 1980) with 100 000 
iterations and restored. The resulting image was then corrected 
for the mean primary beam response of the ATCA antennas. 
Also the separate total-intensity images in both frequency bands 
were formed. Thes e images w ere de-convolved using also the 
CLEAN algorithm ( Clark 1 980) with 1 00 000 iterations and then 
restored and primary beam corrected as the image formed using 
the data in both frequencies before. 

The radio morphology in both frequency bands shows two 
radio arcs (see figure |8]l. There are also two radio bright knots 
within the arcs which do neither have a counterpart in X-rays 
(see figure|9]l nor in the optical (dss2red). 

The flux densities of the sources were obtained by integrat- 
ing the emission within multiple polygons enclosing the rem- 
nant. The rms noise in the image measured by similarly inte- 
grating over nearby source-free regions. The spectral indices 
were computed using the total flux density measurements in 
both available frequencies 1.384 GHz and 2.368 GHz. Fitting 
the data with a power law S oc v", we obtained a spectral in- 
dex for the northern knot a = -0.05 + 0.34 and for the south- 
em a = -0.12 + 0.34, respectively. The radio emission of the 
two arcs, however, h ave a spectral index a = -0.71 + 0.34. 
dPe Hortaet al.ll2012h 

4. Observations in other wavelength regimes 

The region around G308.4-1 .4 was observed in the four infrared 
bands 3, 5, 12 and 22 fim during the Wide-field Infrared Survey 




13''43"' --12"' 'll'" -10"' 

Right Ascension (JKOOO) 



Fig. 8. Grey-scale image of the SNR candidate G308.4-1.4 at 
1.3 GHz. The gray-scale range is -2 to 24 mJy beam"' . 
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Fig. 9. Chandra 20' x 20' image of G308.4-1 .4 overlaid with the 
ATCA data at 1.4 GHz. The contour levels are 2, 3, 4, 6, 8, 15 
mJy beam The radio image has a FWHM of ~ 43" and the 
X-ray image is binned with 1" per pixel. 



(WISE, IWright et al .ll2010h. With the 22 ^im filter an extended 
infrared emission, overlapping with the radio and X-ray emis- 
sion regions, is detectable (see figure [TOli. 

The shock wave of a SN which is running into the interstellar 
medium can interact wi th adjacent molecu lar clouds and produce 
y-rays via decays of jiP dAbdo et al.l201 ih . We therefore checked 
the archive of the Fermi Large- Area Telescope for a y-ray source 
at the position of G308.4-1 .4. The Fermi y-ray satellite monitors 
the sky continuously since August 2008. However, no source 
has been detecte d close to G308.4-1 .4 in the Fermi-L AT Second 
Source Catalog dThe Fermi-LAT Collaborationl201 1() and inves- 
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jgh^jmgQS 41'"00=' 40"30' 

Right Ascension (J2000) 



Fig. 10. 10' X 10' WISE image around G308.4-1.4 at 22 micron 
with the radio contour lines in black at 2, 3, 4, 6, 8 and 15 mJy 
and the X-ray contours in white. 



tigating the images ourselves we found no y-ray emission match- 
ing the X-ray morphology. 

5. Summary and Discussion 

The high resolution imaging and spectroscopic capability of 
Chandra revealed the X-ray emission of G308.4-1.4 as an ex- 
tended plasma. The emission region is spherically symmetric 
with a diameter of 8.4 arc-minutes and an excess emission to- 
wards the east. 

Comparing the fitted hh with the integrated hydrogen 
column density towards the source of Kalberla et al. (2005) 
and iDickev & LockmanI (il990l) places G308.4-1.4 within our 
Galaxy. 

The plasma model which is used to fit the energy distribution 
of G308.4-1.4 was especially designed to improve the modeling 
of X-ray spectra of SNR. Other plasma models assume ioniza- 
tion equilibrium and are therefore not able to reproduce the spec- 
tra of young SNRs which are not ol d enough to reach ionization 
equilibrium (Bork owski et al.ir2001h . Furthermore, the tempera- 
ture of the plasma of approximately 6 million Kevin is in the 
typical range for SNRs (McKee 1987). 

The analysis of the archival ATCA radio data revealed non- 
thermal emission with a spectral index of ^ ~0-7, w hich is typ- 
ical for young to middle-aged SNRs ( lDickel|[l991h . The radio 
morphology at 0.843, 1.432 and 2.448 GHz is characterized by 
two arcs. The eastern arc matches well with the X-ray and in- 
frared contours. The second arc which is bending in the west 
and decreasing brightness towards the center of G308.4-1.4 has 
no counterpart in other wavelengths. Its origin will be discussed 
in more detail in section 153] 

Infrared radiation from SNRs is expected to come either 
from thermal emission of dust grains in the hot plasma which 
are heated through collision or from the forbidden lines of ele- 
ments like Neon, Oxygen and Iron (Seok et al. 2008). Thus, the 
infrared emission matching the radio and X-ray morphology is 
another hint that G308.4-1.4 is indeed a SNR. 



The lack of y-ray emission from G308.4-1.4 can be ex- 
plained by the poor spatial resolution of the Fermi-LAT tele- 
scope (~ 1° at 1 GeV), which makes it hard to distinguish faint 
extended sources from the diffuse emission of the Galactic plane. 
Hence, from the 2 74 known SNRs only 7 have been det ected by 
Fermi-LAT so far (iThe Fermi-LAT Collaborationl201 lb . 

There are various physical scenarios to cause a plasma to 
emit X-rays. However, the two spectral models, APEC and 
MEKAL, which are normally used to describe the X-ray emis- 
sion of Galaxy Clusters (ClGs) cannot reproduce the spectrum of 
G308.4-1.4. Moreover, the temperatur e of ClGs is normally sev- 
eral 10^ K (Bohringer & Wernedl20ll) , one order of magnitude 
higher than the fitted value. Thus, we can conclude that G308.4- 
1.4 is not a CIG. In addition, the possibility of G308.4-1.4being 
a planetary nebula (PN) can be ruled out as the X-ray flux of 
G308.4-1.4 is three orders of magnitude higher than the flux of 
the brightest PN detected in X-rays, which is at a distance of 
» 1 kpc (Guerrero et al..i.2006) . Furthermore, the temperature of 
G308.4-1.4 is a factor of two higher than max imal temperature 
of the hotter, type 2 PNe (iGuerrero et al.ll2000l) . 

Based on this results we conclude that G308.4-1.4 is indeed 
a supernova remnant with its shock wave expanding in more 
dense ISM in the east. Using the inferred spectral parameters we 
can estimate fundamental characteristics of the remnant, such as 
the distance from earth and its age. 

5.1. Distance 

To derive a rough estimate of the distance between SNR 
G308.4-1.4 and Earth we use a method described by 
various authors ("Sun et alj 119991: ISafi-Harb et aP 1200 li: 
iGonzalez & Safi-H arb 2003?): From the dist ribution of m ean 
color excess (Eb-v) per kiloparsec derived by 'Lucke' ('1978') we 
found (Es-v) = (0.3 +0.1) mag kpc"'. Using the fitted hydrogen 
column density for the whole remnant riu = (1.02 + 0.04) x 10^^ 
cm"^, the relation between the hy drogen column den s ity wh 
and the visual extinction Ay of iPredehl & SchmittI (Il995h 
Hfj = (1.79 + 0.03) X 10-' cm^^Ay and adopting the relation 
between Ay and the color excess (-Eb-v) = (3.2 + 0.2)Av 
(Zombeck 2007) the distance to G308.4-1.4 is d ^ 5.9 + 2.0 
kpc. We note, that this distance is probably an underestimation 
of the real distance to G308.4-1.4. This is due to the fact that 
the mean color excess (Eb-v) which we used was derived for 
reddening layer up to 2 kpc. In the case of distance estimates 
for Galactic PNe the reddening inferred distances are mostly 
significantly smaller than distances derived with other methods 
(Phillips 2006). Thus, c/ w 6 kpc is a lower limit for the distance 
of SNR G308.4-1.4. 

Another approach to determine the distance to G308.4- 
1 .4 uses the fact that the radius of the remnant can be de- 
termined in two ways, where depend on d with different ex- 
ponents. The first equation is a result of the Sedov analysis dis- 
cussed in the next chapter (equation|4|i and depends on t/'^^. The 
second equation to determine d is simply the comparison of ra- 
dius of the remnant in pc (Rs) and in arcmin {9), d - R^/O, which 
depends linearly on the distance unit d. Using this dependencies 
we derive 

. = 7420. (. [arcmin])-/^ (^) (^) =9.8-kpc, 

(2) 

where is the post-shock temperature and £51 the SN explosion 
released kinetic energy in units of 10^' ergs. In the following all 
important quantities are given in units of t/g g = ii/9.8 kpc. Using 
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the flux values deduced for G308.4-1.4 we compute its X-ray Table 4. Fundamental parameters of SNR G308.4-1.4. 
luminosity to be Li 



0.5- 
'X 



1. 



J +0.3 
-0.2 



X 10^^ d-' erg/s. 



5.2. Age of remnant 

First we have to deduce the post-shock hydrogen density A^h- 
Using the normalization constant derived in the vgnei-model fit 
and assuming that the hydrogen and electron densities are con- 
stant over the volume V of the emission region we derive the 

following relation: A^eA^H — Norm ■ -j^Ay cm"^. V can be ap- 
proximated by a sphere with radius r w ■ D of which only 
a fraction of ~ is bright enough to be detected in the X- 



1.14 X 10^4 (r[arcmin])3f;3^^cm3 



is the distance to G308.4-1.4 in 



ray regime. Therefore, V 

7.98 X lQ^^dl^cm\ Here dtpc 
kpc. For a fully ionized plasma with cosmic abundances 10% 
He) the ratio between hydrogen and electron density can be es- 
timated by Nh/Ns w 0.8. Thus, the post-shock hydrogen density 

and the corresponding swept-up 



A^H = (0.77 + 0.11) d/gg^^ cm-3 



ll.ldl'^Me (ISlane et al.ll200l . 



mass M = LANumuV 

The best fit to the energy spectrum suggest that G308.4- 
1.4 is in the Sedov-phase. Additionally, the age for the transition 
from the Sedov-pha se into the radiative phase is about 29 kyrs 
dBlondin et al.lll998). Therefore , the age of the remnant can be 
approximated bv (lMcKeelll987h 



f = 2.71 X lO^'l^l ^4„yrs 



(3) 



Assuming that the explosion energy is equal to the canonical 
value of 10^' ergs and the fitted temperature is approximately 
Tfi, the derived value for the hydrogen density suggest that SNR 
G308.4-l.4is * 5000 - 7500 years old. Using this resuh we can 
calculate the theoretical radius of the remnant CCulhaneiil977i) 



R, = 0.34 



f 5 pC 



ii.9!|:^4,pc 



and the shock velocity ( ISedov|[T959l) 



^^= 730!- km/s. 



(4) 



(5) 



which is independent on the assumed distance. 

The Sedov analysis lead to the conclusion that the SNR is 
relatively young and is expanding with a velocity in the order of 
X 730 km/s. All fundamental parameters of the SNR are summa- 
rized in table |4] However, the numbers in this table should only 
be seen as rough estimates as we could not determine whether 
a full temperature equilibrium is akeady established. If not, a 
generally higher ion temperature which determines the SNR dy- 
namics, would lead to a higher shock velocity and thus to an 
e.g. overestimation of the remnant's age. 

5.3. The central point sources CXOU 
J1 341 24.22-634352.0 and CXOU 
J1 341 27.1 2-634327.7 

The X-ray point source CXOU J134124.22-634352.0, located 
close to the geometrical center of the SNR G308.4-1.4, is seen 
to exhibit variable and flaring X-ray emission. It lacks a radio 
counterpart, has no compact/extended X-ray emission as seen in 
young and powerful pulsars and its best-fit spectral model con- 
sists of two blackbodies. Thus, source # 1 shows some proper- 
ties which have been seen in central compact objects (CCO, see 



SNR 


distance 


9.8!^;« kpc 


age 


5000 - 7500 yrs 


radius 


~4.2711.9!]X8 pc 


expansion velocity 


730!|^S km/s 



Got thelf & Halpernll2008l for a review). Moreover, one of these 
CCDs shows flares, IE 161348-5055 in RCW 103, where its 
flux changes between 0.8 and 60 x 10"'^ ergs cm"^ s~'. This is 
comparable to source #1 where the maximal flux is x 10"" ergs 
cm"^ s ' and the upper limit derived from the SWIFT observa- 
tions is w 5 X lO"'-' ergs cm"^ s"'. 

Using the distance derived in chapter ISTI the luminosity of 



CXOU J134124.22-634352.0 is L^^-io 



= 2.5-^:^ X 10^4 d,l 



erg/s. This is s lightly higher than the normal population of CCO 
(lBeckedl2009l p. 121), only the flaring CCO in RCW 103 ha s a 
comparable luminosity of 0.1 - 8 x 10^"* erg/s jde Lucall2008h . 

Additionally, the X-ray-to-optical flux ratio and the fitted «h, 
which is in the order of what has been observed for G308.4-1.4, 
suggest a CCO interpretation for source # 1 . 

Assuming that the source is indeed the compact remnant 
of SNR G308.4-1.4 and using the inferred center and the age 
of the SNR we can derive the proper motion of the object 

ju = yJ(RA, - /?A#i)2 + (DECc - DEC#i)^/t ~ 45 .4" /f. Using 
as minimal age the derived lower limit for f we deduce an upper 
limit for the proper motion of // < 17 mas/yr » 470 km/s. This 
limit is higher than the mean two -dimensional proper motion of 
pulsars in SNRs of ~ 227 km/s (iHobbs et al.l l2005). but would 
not be an exception. For example PSR B1830-08 in SNR W41 
has a 2d proper motion of a: 730 km/s. 

Nevertheless, the lower limit of the inferred emitting radius 
of the blackbody in the BB+BB and PL+BB spectral fit is in the 
order of the expected value for a neutron star of 10 to 20 km 
only if the source is at a distance larger than 14.5 kpc (see fig- 
ure|7]). Moreover, we found an infrared and optical source which 
position is consistent with source # 1 . The spectral distribution 
of the source allows us to fit an APEC model, which ^ten et al.l 
(2005) used to model the energy distribution in the X-ray band 
of the M dwarf flare star EV Lacertae. In addition, the fitted hy- 
drogen column density is consistent to that of a nearby star 

Using the stellar density in direction of G308.4-1.4 we can 
derive the chance association of the observed 2MASS source 
with the compact remnant candidate. The chance association is 
given by f com - J^^J^^ SDEC, where A^ is the number of 
sources detected within a rectangular region of length Ira and 
width loed. and 6RA, SDEC are the errors in the position of the 
source. The 2MASS catalog contains 2206 point-like sources in 
a box with side length of 5 arcmin around the center of the rem- 
nant and the chance association is 1.3%. Thus, the possibility 
of a false association between the X-ray and the 2MASS source 
cannot be excluded. 

Until now, we did not discuss the western radio arc which has 
no counterpart in the X-ray and infrared image. The arc could 
be interpreted as a relativistic radio jet of a source at the geo- 
metrical center of the SNR. If so, the central source is either a 
black ho le of stellar mass or a neutron star in a binary system 
dMirabel & Rodriguez 1999, and reference therein). Especially 
the possibility that the central source is a black hole is an inter- 
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esting speculation which, if proven, would mean that G308.4- 
1.4 is the remnant of a Type II core-collapse SN, and the first 
one known which left over a black hole. 

To obtain a rough estimate of the source flux of the other 
source in the center of the remnant, source # 10, we assumed 
that the source is a CCO, should have at least a spectrum with 
a blackbody of temperature 2.6 million Kevin like the CCO in 
Puppis-A (cf. table 6.4, Becker 2009). Using the WebPIMMSQ 
tool with the source count rate and the fitted nn of the SNR the 
flux in the 0.3 - 3.5 keV range is fx ~ 2.2 x 10"'"* ergs cm^^s"' 
and in the 0.5 - 10 keV range is ~ 1-9 x 10"'"* ergs cm^^s"'. 
The normalization of a blackbody with this flux, temperature and 
rill is 2.5 X 10"^ which corresponds to an emitting radius of the 
source of ^ 9t/g ^ km. This value is in perfect agreement with the 
expected value for the neutron star radius. The corresponding 
luminosity in the energy band 0.5 - 10 keV is L^-^""' = 2.1 x 
10^^ t/g g erg/s, a typical value for CCOs. Furthermore, the X-ray- 
to-visual flux ratio is \og{ fxl fv) > 0.1 1 taking into account that 
no optical source was detected in the USNO- Bl.O catalog an d 
the limiting magnitude of this catalog is 21 (Monet et al.ll2003l) . 
Thus, no clear evidence for the origin of the emission of source 
#10 could be found. 

Finally we mention that the detector support structure of 
Chandra covers the central part of G308.4-1.4 (see figure|2]l. It 
is therefore not excluded from our observations that other faint 
X-ray point sources are located at this position. However, the 
non-detection of a source in the SWIFT data at the assumed po- 
sition of the remnant's expansions center sets a strict limit on the 
flux of such an object. 
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